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Abstract 

We describe the features of the stereoscopic approach in the imaging atmospheric 

Cherenkov technique, and discuss the performance of future low threshold telescope arrays. 
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J!v ! 1. Introduction 

The high detection rate, the ability of effective separation of electromagnetic and hadronic 
showers, and good accuracy of reconstruction of air shower parameters are three remarkable 
features of the imaging atmospheric Cherenkov telescope (IACT) technique (see e.g. Cawley and 
Weekes 1996, Aharonian and Akerlof 1997). The recent reports about detection of TeV 7-rays 
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from several objects by nine groups running imaging telescopes both in the northern (Whipple, 
HEGRA, CAT, Telescope Array, CrAO, SHALON, TACTIC) and southern (CANGAROO, 
Durham) hemispheres confirm the early expectations concerning the potential of the technique 
(Weekes and Turver 1977, Stepanian 1983, Hillas 1985), and provide a solid basis for the 
ground-based gamma ray astronomy (Weekes et al. 1997a). 

In particular, the Whipple 10 m diameter imaging telescope, today's most sensitive single 
telescope, provides detection of more than 100 7-rays from the Crab with ~ 7 standard deviation 
statistical significance within only 1 h observation of the source. Thus the 100 h source exposure 
by this instrument should be enough to reveal point 7-ray sources above 250 GeV at the flux 
level of 0.07 Crab. On the other hand, the '7-sigma-per-l hour' signal from the Crab implies that 
any short outburst of a TeV source at the level exceeding the Crab flux can be discovered during 
the observation time of less than 30 minutes. This important feature of the IACT technique 
was convincingly demonstrated by the Whipple group when two dramatic flares from Mrk 421 
were detected in May 1996 (Gaidos et al. 1996). Significant improvement of the performance 
of this telescope is expected in 1998, after installation a new, 541-pixel high resolution camera 
(Lamb et al. 1995a). 

While the cameras with relatively modest pixel size of about 0.25° provide an adequate 
quality of imaging of air showers produced by TeV 7-rays (Plyasheshnikov and Konopelko 1990, 
Zyskin et al. 1994, Aharonian et al. 1995), a smaller pixel size is preferable (a) for lowering 
the energy threshold, and (b) for extending the dynamical range towards the higher energies 
by observations of 7-ray sources at large zenith angles. Both possibilities recently have been 
demonstrated respectively by CAT (Goret et al. 1997) and CANGAROO (Tanimori et al. 
1997) groups. 

Further qualitative improvement of the IACT technique in the next few years is likely to 
proceed in two (complementary) directions: (a) implementation of the so called stereo imaging, 
and (b) reduction of the energy threshold towards the sub- 100 GeV domain. 
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Fig. 1: Density of Cherenkov light photons in 10 GeV and 100 GeV^-ray showers as a function 
of the radial distance from the shower core. 

2. Stereoscopic Approach 



The concept of stereo imaging is based on the simultaneous detection of air showers in 
different projections by at least two telescopes separated at a distance comparable with the ra- 
dius of the Cherenkov light pool, Rq ~ 100 m (see Fig. 1). Compared to single telescopes, the 
stereoscopic approach allows (i) unambiguous and precise reconstruction of the shower param- 
eters on an event-by-event basis, (ii) superior rejection of hadronic showers, and (Hi) effective 
suppression of the background light of different origins - night sky background (N.S.B.), local 
muons, etc (Aharonian et al. 1993, Aharonian 1993, Krennrich and Lamb 1995, Stepanian 
1995, Chadwick et al. 1996, Aharonian et al. 1997). 

The (only) disadvantage of the stereoscopic approach is a significant loss of the detection 
rate due to the overlap of the shower collection areas of individual telescopes located from each 
other at distances < 2Rq. However, the loss in statistics is partially compensated, especially for 
steep spectra of primary 7-rays, by reduction of the energy threshold of the telescopes operating 
in the coincidence mode. 

While the effective collection area of a single telescope is determined by the radius of the 
Cherenkov light pool, the shower integration area of a multi-telescope system is due to the 
total geometrical area of the array which may be gathered from individually triggered groups 
of telescopes - IACT cells (Aharonian et al. 1997). The compromise between two principal 
conditions, namely (a) detection of 7-ray induced showers by several telescopes of the cell and 
(b) minimization of correlation between images in different telescopes, dictate the linear size of 
the cell: L ~ Rq. The design of the IACT cell depends on the specific detection requirements. 
However, since the accuracy of reconstructed shower parameters continues to be improved 
noticeably up to 3 or 4 telescopes in coincidence (Aharonian et al. 1997), the optimum design 
of the cell seems to be a triangular or quadrangular arrangement of IACTs. The concept of the 
IACT cell makes straightforward and rather general the predictions concerning the performance 
of multitelescope arrays. 



2.1. Energy Threshold 
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Fig. 2: Energy threshold of an imaging Cherenkov 
telescope versus the telescope's aperture S'ph.e.- 
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Fig. 3: Differential detection rates of '-y-rays from 
a point source with differential power-law index 2.5 
and integral flux J 7 (> 100 GeV) = 10~ n ph/cm 2 s, 
as well as detection rates of isotropic cosmic ray 
protons and electrons (^xlO -4 ^ for the 4-fold tele- 
scope coincidence. The curves 2 and 4 are cal- 
culated for a camera with pixel size of 0.15°. All 
other curves are calculated for the pixel size 0.25°. 

to-photoelectron' conversion factor, determined 



The effective energy threshold of the 
imaging telescopes, E th , is basicly deter- 
mined by two conditions: (i) the acciden- 
tal trigger rate due to the N.S.B. should 
not exceed the detection rate of 7-rays, and 
(ii) the number of photoelectrons in image 
should be sufficient for proper image anal- 
ysis; typically nj,™ ~ 100. 

Effective suppression of the N.S.B. in 
the IACT technique is provided by the trig- 
ger criterion which requires signals above 
threshold in several adjacent pixels. Since 
the minimum number of pixels in coinci- 
dence is limited by the characteristic im- 
age size, the high resolution cameras with 
pixel size of about 0.1° have a certain ad- 
vantage from the point of view of lowering 
the energy threshold by using higher trig- 
ger multiplicity. 
Meanwhile, the disadvantage of (more eco- 
nomical) cameras with relatively large (e.g. 
0.25°) pixel size can be effectively com- 
pensated by the requirement of simulta- 
neous detection of a shower by > 2 tele- 
scopes with a local 2-pixel coincidence trig- 
ger condition. Due to the flat lateral dis- 
tribution of the Cherenkov radiation of 7- 
ray showers this requirement does not af- 
fect the 7-ray detection efficiency if the 
distance between telescopes does not ex- 
ceed 100 m. Assuming now that the lo- 
cal and the system triggers eliminate ef- 
fectively the non- air- shower backgrounds, 
the threshold E t h could be estimated from 
the condition (ii), namely from the equa- 
tion n^ ] = p(E) x S'ph.e, where p(E) is 
the density of the Cherenkov photons at 
plateau produced by a primary 7-ray of 
energy E, and S'ph.e. is the telescope aper- 
ture for the number of detected photoelec- 
trons, Oph.e. = >->mir ' Xph— »e- Here O m i r 

is the geometrical area of the optical re- 
flector and Xph->c is the so-called 'photon- 
by quantum efficiency of the light detec- 



tors, the reflectivity of the mirror, etc. If for order-of-magnitude estimates one assumes 
that p(E) ~ PoE, the condition ni™™ = 100 gives, for example at 2 km above see level, 
E t h = 100 (S'ph.c./lOm 2 )^ 1 GeV. In Fig. 2 we show accurate calculations of the dependence of 
E th on the telescope aperture S p ^ e , for different altitudes. 

From Fig. 2 one may conclude that the energy threshold around 100 GeV can be achieved, 
in principle, by a system of 10 m diameter IACTs equipped with conventional PMT based 
cameras with typical Xph-^c — 0.15. The Monte Carlo calculations of detection rates for quad- 
rangular 100 m x 100 m cell consisting of four IACTs with aperture S p h. e . = 10 m 2 show that 
the detection rate of 7-ray showers, with the cores inside the cell, indeed peaks at 100 GeV 
(Fig. 3). The efficiency of detection of 100 GeV photons even in the stringent 4-fold telescope 
coincidence mode is about 50%, and reaches the 100% level already at 200 GeV. 



Note that n 



(min) 
ph.c. 



100 is a rather conservative condition. For relatively strong 7-ray 
sources, when a comprehensive image analysis is not required, the energy threshold can be 
lowered down to 50 GeV. 



2.2. Reconstruction of Shower Parameters 



Cherenkov light 
image 



The procedure of reconstruction of the arrival direction and core position of 7-ray showers 
by 4 telescopes of the cell is demonstrated in Fig. 4. Although one pair of detected images is 
sufficient to determine unambiguously the shower parameters (except for the events close to 
the line connecting two telescopes), the information from other images significantly improves 
the accuracy of estimated shower parameters. 

The one dimensional angular distributions 
of showers produced by 7-rays from a point 
source with energy E — 100 GeV, 300 GeV, and 
1 TeV are shown in Fig. 5a. The two-dimensional 
distribution of the reconstructed directions of 7- 
rays from a point source with E~ 2 spectrum is 
shown in Fig. 6. The values of the 'la' angu- 
lar resolution (68% acceptance of 7-ray showers), 
are presented in Table 1. 

Detection of the shower images by > 2 
telescopes allows precise determination of the 
shower core position (impact parameter). The 
distributions of uncertainties in reconstructed 
impact parameter for different primary energies 
are shown in Fig. 5b. The half width at half 
maximum (HWHM) values derived from these 
distributions are presented in Table 1. 

Determination of the shower core position 
with accuracy < 10 m is quite sufficient, due to 
the flat lateral distribution of the Cherenkov light from the 7-ray showers (see Fig. 1), for the 
accurate evaluation of the energy using a relation between the primary energy E and the total 
number of photoelectrons in the image. The accuracy of the energy estimate is essentially 
due to the fluctuations which are large (> 30%) for very close (R < 60 m) as well as distant 
(R > 120 m) showers (see Konopelko 1997). However, most of the showers required to be 
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Fig. 4-' Reconstruction of the shower param- 
eters by the J^-IACT cell. 
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Fig. 5: a (left): One- dimensional angular distribution of'y-rau showers, b (right): Uncertainties 
in the reconstructed impact parameter of 7 -ray showers. 

detected by all four telescopes are located between 60 m and 120 m from at least 2 telescopes 
of the cell. Therefore the fluctuations appear to be small enough to provide determination of 
the energy of 7-ray showers with an accuracy of about 20% in the broad energy range from 
100 GeV to 1 TeV (see Table 1). 

2.3. Gamma/Hadron Separation 
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Determination of the arrival direction 
of individual 7-ray showers with an accu- 
racy < 0.1° allows effective suppression of 
the cosmic ray background for point sources 
by a factor of k^ = (259 /P si) 2 , where \1> 
is the the angular diameter of the 'trigger 
zone'. For point sources typically \I/ ~ 3°, 
thus /4r ir ' ) — 1/300. Further improvement 
of the signal-to-noise ratio is provided by ex- 
ploitation of the intrinsic differences between 
the development of the electromagnetic and 
hadronic cascades. These differences, caused 
by the large transverse momenta of the sec- 
ondary hadrons, the deeper penetration of the 
hadronic cascades into the atmosphere, and 
essential fluctuations in the hadronic cascades 
(e.g. Hillas 1996), are effectively transferred 
to to the shape of the Cherenkov light im- 
ages. The 7-ray images have compact shape 
and regular structure, in contrast to rather 
diffuse and chaotic images of the hadronic showers. For single telescopes this provides rejection 
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Fig. 6: Two dimensional distribution of re- 
constructed arrival directions of 'y-rays from 
a point source with the differential spectrum 
dJ^/dE ex E~ 2 . Circles correspond to the 
30%, 50%, 70%, 85% •y-rays acceptance. 
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Fig. 7: Cherenkov light images of vertical 100-GeV / ~f-ray (left) and 300-GeV proton (right) 
detected by four telescopes of a 100 m x 100 m I ACT cell. Both showers cross the center of the 
cell. While the images of the proton shower in two telescopes look like ^-ray showers (narrow 
and regular), the images obtained by two other telescopes clearly indicate the hadronic nature 
of the shower. 

of hadronic showers, based on the so-called shape cuts, by a factor of 10 or so (see e.g. Fegan 
1997). 

Determination of the arrival direction, energy and core location of showers on an event- 
by- event basis with accuracy < 20%, < 0.1°, and < 10 m, respectively, allows an optimization 
of the image 'cuts' which separate the CR and 7-ray domains in the shower parameter space, 
and thus improves the rejection efficiency for hadronic showers. Furthermore, even though 
the images in different projections are not entirely independent of each other, the correlation is 
only partial. Therefore, the joint analysis dramatically improves the gamma/hadron separation 
power (this effect is demonstrated in Fig. 7). The characteristic values of the acceptance of 
7-ray (k 7 ) and proton {Kh} induced showers, as well as the so-called Q-factor, Q = ft 7 /% (the 
enhancement of the signal-to- noise ratio after the application of the cuts; see, e.g. Fegan 1997), 
are presented in Table 2. 

Table 1: The angular, distance, and energy resolutions of the J^-1ACT cell. 
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It is seen that rejection efficiency of hadronic showers significantly increases with the 
number of detected images, and under the condition of 50% acceptance of 7-ray events the 
rejection factor k^ 1 of 4-IACT cell can be as large as 300. This results in the Q-factor of about 
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Fig. 8: Detection rates of '-/-rays (solid lines), as well as cosmic ray protons (dashed) and 
electrons (dot-dashed) by the 100 m x 100 m IACT cell, after application of software image 
cuts, (a) - point source, and (b) - '1° ' diffuse source of / y-rays. In both cases an integral 
flux of t -rays J 7 (> 100 GeV) = 10~ n ph/cm 2 s with power-law indices a = 2 and a = 3 is 
assumed. In the case of point sources the combined 50% acceptance for '-/-rays is achieved by 
applying directional and relaxed shape cuts, while in the case of a 1° ' source the same 50% 
7 -ray acceptance is a result of tight shape cut. 

8. Note that such effective gamma/hadron separation is provided by application of a 'scaled' 
Width parameter (Konopelko 1995). The other parameters, like image Length, time structure, 
UV-content of the Cherenkov light flashes, etc., are less effective, but they still can improve the 
efficiency of suppression of the hadronic showers being applied separately. However, in the case 
of stereoscopic observations, the Width cuts in different projections are so effective, that they 
practically do not leave a room for further enhancement of the signal (Aharonian et al. 1997). 
The suppression of the cosmic ray background based on the directional and shape criteria 



for point sources could be as effective as k ct = k 



(dir) 



K h 



10 with the Q-factor close to 100. For 



comparison, CR rejection efficiency and Q-factor in case of single IACTs, after the application 
of both orientational (Alpha) and shape cuts, hardly can exceed 10~ 3 and 10, respectively. 

Due to the high rejection efficiency of hadronic showers by the IACT cell, the cosmic ray 



Table 2: CR background rejection efficiency of IACT 'cell' using parameter Width with and 
without 'scaling'. In the analysis the N-fold Width cuts were used (N = 2,3,4; N is the number 
of the images passed the cut). 



Parameter: 


Scaling: 




2 


3 


4 


Width 


No 


Q-factor 


0.792 
0.116 
2.32 


0.585 
0.038 
3.00 


0.489 
0.016 
3.87 


Width 


Yes 


K/sy 

K h 

Q-factor 


0.891 
0.076 
3.23 


0.634 
0.017 
4.86 


0.490 

0.0036 

8.33 



background at energies below 200 GeV is dominated by showers produced by electrons (see 
Fig. 8 a,b). This allows detection of 7-rays from point sources under almost background-free 
conditions if the photon flux above 100 GeV exceeds 10 -11 ph/cm 2 s (approximately 1/20 of the 
Crab flux expected at these energies). 

3. HEGRA: The Stereo Imaging Does Work ! 
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The recent observations of the Crab 
Nebula and Mrk 501 by HEGRA stereoscopic 
system of 4 IACTs (Daum 1997, Herman 
1997, Hofmann 1997) confirm the early pre- 
dictions for the performance of the instrument 
(Aharonian 1993). Fig. 9 shows the angular 
distribution of showers detected from the di- 
rection of the Crab, and selected after the 
image shape cuts (Daum et al. 1997). The 
confinement of the 7-ray signal from a point 
source in a very small angular region (with ra- 
dius ~ 0.1°) of the available two-dimensional 
phase space, in addition to the significant sup- 
pression of hadronic showers at the trigger 
level, results in a strong 4<r-per-l h signal al- 
ready before the gamma/hadron separation. 
The shape cuts provide further, by a factor of 
100, suppression of the hadronic background, 
while maintaining the efficiency of the accep- 
tance of 7-ray events at the level of 50%. This 



Fig. 9: Distribution of the reconstructed direc- 
tions of the showers detected from the Crab by 
the HEGRA I ACT system (from Daum et al. 
1997). 

leads to only ~ 1 background event in the 'signal' region while the rate of 7-rays from the Crab 
exceeds 20 events/h. This enables (i) an effective search for VHE 7-ray point sources with fluxes 
down to 0.1 Crab at almost background free conditions, and thus allows drastic reduction of the 
observation time (by a factor of 10) compared to four independently operating telescopes, (ii) 
detailed spectroscopy of strong (Crab-like) VHE emitters, (in) study of the spatial distribution 
of 7-ray production regions on arcminute scales, (iv) effective search for extended sources with 
angular size up to ~ 1° at the flux level of 0.1 Crab. The exploitation of the 'nominal' 5- 
IACT HEGRA system with improved trigger condition should allow 5a detection of faint 7-ray 
sources at the flux level of 0.025 Crab (Aharonian 1997, Hofmann 1997) which corresponds to 
~ 10~ 12 erg/cm 2 s energy flux above the effective energy threshold of the instrument of about 
500 GeV. 



4. Future 100 GeV Class IACT Arrays 



One of the principal issues for future detectors is the choice of the energy domain. If one 
limits the energy region to a relatively modest threshold around 100 GeV, the performance 
of IACT arrays and their practical implementation can be prediced with high confidence. In 



practice, an energy threshold of 100 GeV can be achieved by a stereoscopic system of IACTs 
consisting of 10 m diameter optical reflectors equipped with conventional PMT-based high 
resolution cameras. Two such arrays of 10 m class telescopes - VERITAS (Weekes et al. 
1997b) and HESS (Hofmann 1997) - are likely to be constructed in the foreseeable future. 



FUTURE GROUND-BASED GAMMA RAY DETECTORS 
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operation modes: 



(1) at small zenith angles 

3 independent stereoscopic systems 
of IACTs (at E< 100 GeV) 

(2) at large zenith angles 

individual IACTs from different 
cells - in coincidence (at E > 100 GeV) 



Fig. 10: Two possible arrangements of (a) 100 GeV class and (b) 10 GeV class IACT arrays. 



After reaching the maximum possible suppression of the cosmic ray background by si- 
multaneous detection of the air showers in different projections - limited basically by intrinsic 
fluctuations in cascade development - the further improvement of the flux sensitivities for a 
given energy threshold can be achieved by increasing the shower collection area, i.e. the num- 
ber of IACT cells. A possible arrangement of a cell-structured array is schematically shown in 
Fig. 10a. For the arrays consisting of a large number of cells, e.g. uq > 10, and for a given 
observation time to, the minimum detectable flux of 7-rays from a source of angular size <fi is 
determined essentially by the product A = no x to from the following two conditions: (1) high 
confidence level of a 7-ray signal, i.e. a = R^(R P + R c )~ 1 / 2 A 1 I 2 G~ 1 (0) > a , and (2) adequate 
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Fig. 11: Minimum detectable fluxes of y -rays from a point source (solid lines) and 1° source 
(dashed lines) within 100 h observation time by an array consisting of twenty five 10 m diameter 
I ACT cells. The calculations correspond to the requirement of a 5a y-ray signal with at least 25 
detected photons. For comparison, we show the following fluxes: (i) power-law extrapolations 
of a 'standard' EGRET source with integral indices (a — 1) = 1 and 2.5; (i) the level of the 
measured ■y-ray fluxes of the Crab Nebula above 300 GeV, assumed in the form of Joab(> 
E) ~ 1.5 x 10~ n (E/l TeV)~ L5 ph/cm 2 s. The minimum detectable fluxes ofy-rays by the single 
HEGRA and Whipple telescopes (indicated as H(l) and W, respectively), and by the HEGRA 
IACT system (H(S)), calculated for 100 h observations of a point source, are also shown. 

photon statistics, i.e. iV 7 = R^A > N . Here i? 7 , _R P , and R e are the detection rates of 7-rays 
from a point source, and the detection rates of cosmic ray protons and electrons by one IACT 
cell, respectively (see Fig. 8a). The parameter G(0) ~ max[l,0/<5#] takes into account the 
angular size of the source. In Fig. 11 we show the minimum detectable 7-ray fluxes from a 
point source and 1° source requiring > 5a appearance of a signal, provided that the number of 
detected 7-rays exceeds 25. Note that at a large exposure, e.g. A > 2500 h, which can be real- 
ized, for example, by 100 h observations of a source with an array consisting of 25 IACT-ce//s, 
</ 7 nm (> 100 GeV) ~ 10 _12 ph/cm 2 s, or the corresponding energy flux /□ ~ 10~ 13 erg/cm 2 s. 

At a small exposure, e.g. A = uq x t < 100 h, the minimum detectable fluxes are limited, 
especially at high energies, by low 7-ray statistics. This could be treated as a disadvantage of 
cell-structured arrays with fixed linear size of the cell, since such arrays with L ~ Re ~ 100 m 
are optimized, to a large extent, to the energy region close to the threshold. Indeed, although 
the homogeneous multi- cell arrays provides rather economical coverage of large detection areas 
at the threshold energies (each telescope is exploited by four cells), at high energies this concept 
becomes less effective. Namely, it cannot use the advantage that the showers initiated by high 
energy photons, E ^> E th , can be detected at distances well beyond 100 m. However this 
disadvantage of homogeneous cell-structured arrays can be compensated, at least partially, by 
observations at large zenith angles by > 2 telescopes in coincidence from different peripheries 



of the array. 
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Fig. 12: Distributions of the reconstructed arrival 
directions of '7 -rays from sources with different an- 
gular radius. The histogram 1 corresponds to the 
two-dimensional distribution shown in Fig. 6. 



From this point of view it would be 
important to equip some of the telescopes 
of the array, located at large distances from 
each other, by very high resolution cameras 
with pixel size ~ 0.1°. 

The determination of arrival direc- 
tions of 7-rays by IACT arrays on an event- 
by-event basis with accuracy < 0.1° makes 
possible a study of the spatial distribution 
of VHE 7-ray sources in arcminute scales. 
This is demonstrated in Fig. 12, where the 
distributions of the angular distance of re- 
constructed directions of 7-rays from the 
center of extended source are shown. From 
this figure we may conclude that for suffi- 
ciently high photon statistics (say > 10 3 
detected 7-rays) it would be possible to 
measure the size of an extended source of 
< 0.1° with accuracy < 1 arcmin. 
20% enables reasonable spectrometry with 



Good energy resolution of about SE/E 
IACT arrays with ability, for example, to reveal spectral features like sharp cutoffs in the 
spectra of VHE sources (see Fig. 13). 

Due to very effective rejection of the hadronic showers by the stereoscopic IACT systems, 
the cosmic ray background at energies close to 100 GeV is dominated by the showers produced 
by primary electrons. This component of electromagnetic showers remains as a part of the 
background which principally cannot be removed, and thus it is the most serious limiting factor 
of flux sensitivities, especially for extended sources. On the other hand, these electromagnetic 
showers with known flux and spectrum of the primary electrons, measured up to energies 1 
TeV (Nishimura et al. 1980), can be used for absolute energy calibration of the IACT arrays. 
Although the current uncertainties in 100 GeV electron fluxes could be as large as factor of 2, 
the future balloon and satellite measurements (e.g. Torii et al. 1997, Miiller et al. 1997, Ormes 
et al. 1997), should be able to reduce significantly the flux uncertainties. 

The energy calibration by using the cosmic ray electrons provides unique tool for the 
continuous (on-line) control of the characteristics of IACT arrays (e.g. the energy threshold, 
the detection area, etc.) during the 7-ray observations. For example, 10 min exposure will be 
enough for detection of more than 100 electrons within the field of view of 100 GeV threshold 
array consisting of 9 IACT cells (see Fig. 8b). It is difficult to overestimate the significance of 
such calibration and control, especially for the study of the spectral characteristics of highly 
variable 7-ray sources on sub-hour timescales. 



5. Sub-100 GeV Ground-Based Detectors 



The strong scientific motivations to fill the gap between the space-based and ground based 
observations recently stimulated several interesting proposals for extending the atmospheric 
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Cherenkov technique down to 10 GeV region (see e.g. Lamb et al. 1995b). Among a variety of 
competing projects this goal could be addressed by the stereoscopic systems of IACTs with the 
telescope apertures of S p h.e. > 50 m 2 . With development of the technology of novel, fast optical 
detectors of high quantum efficiency (Xph-*c > 50%) and design of (relatively inexpensive) 20 m 
class reflectors with an adequate optical quality for the Cherenkov imaging on scales of several 
arcminutes, it would be possible to reduce the detection threshold down to 20 GeV, or even 10 
GeV for systems installed at high mountain elevations (e.g. 3.5 km a.s.l). Considerable R&D 
efforts, in particular by the Munich group (Lorenz 1997, Mirzoyan 1997) are already started 
in both directions. It is expected that this activity will result in construction of a large single 
reflector telescope (MAGIC) which, together with two ongoing projects of low energy threshold 
air Cherenkov detectors based on large mirror assemblies of the existing solar power plants - 
STACEE (Ong et al. 1997) and CELESTE (Pare et al. 1997) - will start to explore the energy 
region below 100 GeV. Moreover, if successful, the MAGIC telescope could be considered as 
a prototype for the basic element in future (post 'VERITAS/HESS') arrays of 10 GeV class 
IACTs. 

The efficiency of the imaging tech- 
nique at energies well below 100 GeV seems 
to be lower than at higher energies. In- 
deed, at energies close to 10 GeV, the 7-ray 
images become less elongated (more cir- 
cular) and less regular. In practice this 
means an introduction of significant un- 
certainties in the reconstruction of the ar- 
rival direction of 7-rays, as well as in the 
gamma/hadron separation strategies. In 
particular, our preliminary results show 
that the angular resolution of stereoscopic 
systems at 10 GeV hardly can be better 
than 0.2°. Also, due to large fluctuations 
in low-energy hadronic showers (especially 
due to the higher probability of production 
of a leading 7r° meson which takes almost 
all energy of the primary hadron) , an essen- 
tial fraction of the cosmic rays could "pass" 
the shape cuts, i.e. could be accepted as 
7-rays. On the other hand, since the elec- 
tron/proton ratio of cosmic rays increases 
at lower energies as oc E~ 0M , the weak- 
ened rejection efficiency of hadronic show- 
ers seems to be still sufficient for operation 
of the stereoscopic sub-100 GeV IACT sys- 
tems in the regime when the background 
is dominated by the cosmic ray electrons. 
This makes straightforward the prediction 
for the flux sensitivity of such systems. In- 
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Fig. 13: Evaluation of differential spectra of 7- 
rays detected by an array for A = n x £ = 50 h. 
The primary spectra (solid lines) have a power- 
law form oc E~ 2 - 2 with different values of a cutoff 
energy: 1 - E* > 10 TeV (corresponding recon- 
structed fluxes are shown by filled quadrangles); 2 
- E* = 1 TeV (reconstructed fluxes - filled trian- 
gles); 3 - E* = 300 GeV (reconstructed fluxes - 
filled circles). The CR background events, remain- 
ing after application of image 'cuts ', are shown by 
open circles. 



deed, since the showers caused by electrons are quite similar to the 7-ray showers, the condition 



of detection of 5a signal of 7-rays with differential spectrum dJ 1 /dE at any energy E within 
the interval [E — AE, E + AE] results in 
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where dF e /dEdQ ~ 4-lO~ 2 E~ 3 ' 2 cm _2 s _1 sr _1 GeV _1 is the measured spectrum of the cosmic ray 
electrons; S e g (E) is the effective area of detection of electromagnetic showers, Kg and /t c .- m . are 
the acceptance of electromagnetic showers after application of the directional and shape cuts; 
AQ ~ no- 2 is the solid angle corresponding to the (lcr) angular resolution of the instrument 
cg(E) at energy E; and t is the exposure time. For characteristic values AE = (5E/E)E ~ 
0.25E, K c .- m . = 0.5, Kg = 0.68, we obtain 
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Our preliminary calculations show that with a stereoscopic system consisting of several 
sub-100 GeV threshold IACTs, we may expect a e ~ 0.2°, and S eS ~ 4 x 10 4 m 2 at 20 GeV 
which results in the integral flux sensitivity J 7 (> 20 GeV) ~ 5 x 10~ n ph/cm 2 s. 

Thus, a few sub-100 GeV imaging telescopes being combined in a stereoscopic system 
could effectively intervene the domain of satellite-based gamma ray instruments like GLAST 
(Bloom 1996). But, of course, the scientific goals of these instruments are quite different. If 
the GLAST with its large (almost 27r) FoV could provide very effective sky surveys, sub-100 
GeV ground based instruments have an obvious advantage for the search and study of variable 
sources. On the other side, sub-100 GeV IACT systems could operate also in the TeV domain 
by observing the 7-ray sources at large zenith angles. If high efficiency of this technique, that 
already has been demonstrated by the CANGAROO group at multi-TeV energies, could be 
extrapolated to the 0.1-1 TeV region, the use of an array consisting of few 10 GeV class IACT 
systems (with 3 or 4 telescopes in each) could cover very broad energy range extending from 
10 GeV to 1 TeV. Note however that more detailed Monte Carlo studies are needed for definite 
conclusions about the performance of sub-100 GeV stereoscopic IACT systems. 

A possible arrangement and operation modes of an array of 10 GeV class telescopes are 
described in Fig. 10. 

6. Summary 

The stereoscopic approach in the imaging atmospheric Cherenkov technique provides su- 
perior gamma/hadron separation power (k,^ < 10~ 2 ), excellent angular resolution (<jg < 0.1°), 
good energy resolution (SE/E ~ 20%), and impressive flux sensitivity (~ 10~ 13 erg/cm 2 s) that 
makes the multi telescope arrays as very effective tools for the study of the sky in 7-rays in a 
broad energy region from 10 GeV to > 10 TeV. The recent results obtained by the HEGRA 
system of imaging telescopes, and by the prototype of the Telescope Array (Teshima 1997) 
generally confirm the early predictions concerning the performance of the technique in the TeV 
regime. The forthcoming HESS and VERITAS arrays of 10 m class imaging telescopes, with 
energy threshold of about 100 GeV or less, will improve the flux sensitivities of current in- 
struments by an order of magnitude, and will enable an effective study of distant extragalactic 



sources with redshifts up to z ~ 1. With development of new technologies of construction of 
20 m class optical reflectors and novel high quantum efficiency (> 50%) optical detectors, it 
will be possible to extend the energy domain of the stereoscopic IACT arrays down to 10 GeV. 
However, it should be emphasized that if the performance of the 100 GeV threshold IACT 
arrays, and their practical implementation can be predicted with high degree of confidence, the 
realization of 10 GeV threshold IACT arrays still remains as an exciting challenge. 

We thank W. Hofmann, A.V Plyasheshnikov, and H.J. Volk for many fruitful discussions. 
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